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ABSTRACT 

The joint likelihood of observable cluster signals reflects the astrophysical evolution of the coupled 
baryonic and dark matter components in massive halos, and its knowledge will enhance cosmologi- 
cal parameter constraints in the coming era of large, multi-wavelength cluster surveys. We present 
a computational study of intrinsic covariance in cluster properties using halo populations derived 
from Millennium Gas Simulations (MGS). The MGS are re-simulations of the original 500 h^^ Mpc 
Millennium Simulation performed with gas dynamics under two different physical treatments: shock 
heating driven by gravity only (GD) and a second treatment with cooling and preheating (PH). We 
examine relationships among structural properties and observable X-ray and Sunyaev-Zel'dovich (SZ) 
signals for samples of thousands of halos with M200 > 5 x 10^^ M© and z < 2. While the X-ray 
scaling behavior of PH model halos at low-redshift offers a good match to local clusters, the model 
exhibits non-standard features testable with larger surveys, including weakly running slopes in hot gas 
observable-mass relations and ~ 10% departures from self-similar redshift evolution for 10^^ M© 
halos at redshift z ~ 1. We find that the form of the joint likelihood of signal pairs is generally 
well-described by a multivariate, log-normal distribution, especially in the PH case which exhibits less 
halo substructure than the GO model. At fixed mass and epoch, joint deviations of signal pairs dis- 
play mainly positive correlations, especially the thermal SZ effect paired with either hot gas fraction 
(r = 0.88/0.69 for PH/GD at z = 0) or X-ray temperature (r = 0.62/0.83). The levels of variance 
in X-ray luminosity, temperature and gas mass fraction are sensitive to the physical treatment, but 
offsetting shifts in the latter two measures maintain a fixed 12% scatter in the integrated SZ signal 
under both gas treatments. We discuss halo mass selection by signal pairs, and find a minimum mass 
scatter of 4% in the PH model by combining thermal SZ and gas fraction measurements. 
Subject headings: galaxies: clusters - cosmology:theory 



1. INTRODUCTION 

Accurate cosmology using surveys of clusters of galax- 
ies requires a robust description of the relations between 
observed cluster signals and underlying halo mass. Even 
without strong prior knowledge of the mass-signal rela- 
tion, cluster counts, in combination with other probes, 
add useful constraining power to cosmological parame- 
ters (jCunha et al.l 120091 ). However, significant improve- 
ments c an be realized when the er ror in mass varia nce is 
known f Lima fc Hul [2001 ICunha fc EvraTdl [20091 ) . Im- 
provements can also be gaine d by extendin g the model 
to multiple observed signals (jCunhal l2008f l. especially 
when an underlying physical model can effectively reduce 
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the dimensional ity of the parameter sub-space associated 
with the model ijYounger et al.ll2006D . The coming era of 
multiple observable signals from combined surveys in op- 
tical, sub-mm and X-ray wavebands invites a more holis- 
tic approach to modeling multi-wavelength signatures of 
clusters. 

Signal covariance characterizes survey selection, in 
terms of mass and additi onal observables. F or the case of 
X-ray selected samples, iNord et al] (|2008f ) demonstrate 
that luminosity-temperature covariance can mimic ap- 
parent evolution in the luminosity-mass relation under 
analysis that combines deep. X-ray-flux limited samples 
with local, shallow ones. 

Employing a selection observable with small mass vari- 
ance minimizes such errors. Recent work has shown 
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that the total gas thermal energy, Y, obs ervable via 
an in tegrated Sunyaev-Zeldovich (SZ) effect (jCarlstromI 
|2004|) or via X-ray imaging and spectroscopy, is a sig- 
nal that scal es as a power-law in mass with only ^ 
15% scatter ([W hite et al.' '2002; Kravt sov efall 120061: 
MaughanI 12007^ Q'Hara et al.. ,2007: Zhang et all l2008l : 



Jeltema et al. However, unbiased estimates of the 



mass selection function for Y or any other signal requires 
accurate knowledge of how the signal-mass scaling rela- 
tion evolves with redshift. The redshift behavior of sig- 
nals is generally not well known empirically, although re- 
cent work has begun to probe evolution in X-ray signal s 
to z ~ 1 (|Maughan et al.l 120061 : IVikhlinin et all l2008f ). 
Emerging samples from wide-area SZ surveys should dra- 
matically improve this situation. 

One can address signal-mass covariance using hydro- 
static, virial or lensing mass estimates from observations, 
but several sources of systematic and statistical error 
challenge this approach. For hydrostati c masses, early 
gas dynamic simulations (lEvrardl 119901 : iNavarro et al.l 
[1996; Thomas 1997) suggested that turbulent gas mo- 
tions drove hydrostatic masses to underestimate true val- 
ues by ~ 20%. More sophisticated recent models, with a 
factor thousand improvement in mass resolution, demon- 
strate this effect at a similar level in the mean, with 
^ 15% scatter among individual syst ems (jRasia et al.l 
[20061: lNagail[2006l : iJeltema et al.ll2008[ ). Cluster masses 
can also be measured by the shear induced on back- 
ground galaxies due to gravitational lensing. With this 
method, individual cluster masses have mass uncertain- 
ties ~ 20% due to cosmic web confusion (Hockstra 2003; 
Ide Putter fc Whit3l2005[ ). but large samples can reduce 
the uncertainty in the mean. 

The mean scaling behavior of samples binned in some 
selection signal offers another empirical path to measur- 
ing covariance. Non-zero covariance between the selec- 
tion and an independent, follow-up signal implies that 
the selection-binned scaling relation of the foUowup sig- 
nal with mass need not match that signal's intrinsic mass 
scaling. Comparison of scaling relations from differently 
selected sani p les th ereby offers insight into covariance. 
iRvkoff et al.l ((2008) offer a first attempt at this exer- 
cise for X-ray luminosity and optical rich ness using the 
optic ally-selected SDSS maxbcg sample ([Koester et al.l 
|2007() . The sample contains ~ 13,000 clusters for which 
weak lensing mass estimates have been made by stacking 
the shear of richness-binned s ub-sam ples (Sheldo n et al.l 
120071 iJohnston et al.ll200l . [Rvkoff et al (2008.') stack 
Rosat All-Sky Survey data (|Voges et al.l 11999^ in the 
same maxbcg richness bins, and find that the mean 
X-ray luminosity-mass relation derived with richness 
binning is consistent at the ~ 2 a level with relations 
derived solely from X-r ay data (|Reiprich fc Bohringeil 
I200I: IStanek et al.ll200l . 

A theoretical approach to studying cluster covari- 
ance is to realize populations via numerical simulation. 
While high resolution treatment of astrophysical pro- 
cesses, including star formation, supernova and AGN 
feedback, galactic winds and therma l conduction have 
been included in recent si mulations ("Dolag et al.l 120051: 
Borgani et al. 2006; Kravt sov et al.l 2006; Si iacki et al.l 
20081: IPuchwein et al.ll2008i ). the computational expense 



has limited sample sizes to typically a few dozen ob- 
jects. A detailed study of population covariance requires 



larger sample sizes, as can be generated by lower reso- 
lution simulations of large cosm ological volumes using a 
more limited physics treatment ([B rvan fc Nor manlll998l : 
iHallman et alll2006l: [Cottlober fc Ycpes 2007). 

We take the latter approach in this paper, focusing on 
the bulk properties of massive halos identified the Millen- 
nium Gas Simulations (MGS), a pair of r esimulations of 
the o riginal 500h~^ Mpc Millennium run (|Springel et al.l 
l2005l ). each with 10^ total particles, half representing gas 
and half dark matter. The pair of runs use different treat- 
ments for the gas physics — a gravity-only (GO) simula- 
tion, sometimes called "adiabatic" , in which entropy is 
increased via shocks, and a simulation with cooling and 
preheating (PH). The former ignores galaxies as both a 
sink for baryons and a source of feedback for the hot intr- 
acluster medium (ICM) . The latter also ignores the mass 
fraction contribution of galaxies, but it approximates the 
feedback effects of galaxy formation by a single param- 
eter, an entropy level imposed as a floor at high red- 
shift (lEvrard & Henrv 1991; Kaiser 1991: Bialek et aP 
[200lt IKay et al. 2007; Gottlober fc Yepcs 2007). Our 
study focuses on samples of ~ 5000 halos with mass 
M > 5 X 10^"^ Mq examined at multiple epochs cov- 
ering the redshift range < z < 2. 

The paper is organized as follows: in Section [2] we dis- 
cuss the details of the Millennium Gas Simulations and 
our halo finding approach. In Section[3l we present mean 
scaling relation behavior in both mass and redshift. We 
then turn to covariance about the mean scaling relations 
in Section m Unless otherwise noted, our units of mass 
are lO^^ft-^^M© and halo mass is defined within a sphere 
encompassing a density contrast Ac — 200 times the crit- 
ical density. 

2. SIMULATIONS 
2.1. Millennium Gas Simulations 

The Millennium Gas Simulations (hereafter MGS) 
are a pair of resimu lations of the original Millennium 
(|Springel et al.l l2005f ). a high-resolution, dark-matter- 
only simulation of a 500 Mpc volume. The sim- 
ulations were run with GADGET-2, treating the gas 
dynamics with smoothed particle hydrodynamic s (SPH) 
(|Springe]|[2005[ ). As described in iHartlev et all (|2008i) . 
the MGS use the initial conditions of the Millennium 
simulation, with 5 x 10® dark matter particles, each of 
mass 1.42 x 10^° h''^ Mq, and 5 x 10® SPH gas par- 
ticles, each of mass 3.12 x 10^ Mq, resulting in a 
mass resolution about 20 times coarser than the orig- 
inal Millennium N-body simulation. The gravitational 
softening length is 25 kpc. The cosmological pa- 
rameters match the original: f2„i = 1 ^ = 0.25, 
fib = 0.045, h = 0.73, and erg = 0.9. While some dif- 
ferences between the simulations are expected due to the 
differe nce in mass resolution and gravitational softening 
length, IHartlev et al.l (j2008[ ) verify the positions of dark 
matter halos to within 50 kpc between the original 
Millennium and the MGS. The value of cts is higher than 
the WMAP3 value ([Spergel et al.ll2006t l. but we do not 
expect it to strongly affect the results of this paper. 

In this paper, we consider two models of the MGS: a GO 
simulation where the only source of gas entropy change 
is from shocks, and a PH simulation with preheating and 
cooling along with shock heating. The GO simulation 
is useful as a base model that can be easily compared 
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to previous hydrodynamic simulations of galaxy clus- 
ters. Furthermore, comparing gravity-only simulations 
to observations highlights the cluster properties that are 
strongly affected by astrophysical processes beyond grav- 
itational heating. 

While adiabatic simulations match the self-similar 
prediction, L ~ T^, for the X-ray luminosity- 
temperature relation, obs er vations show a steep e r slop e 
(jArnaud fc EvrardI [l999l : lOsmond k PonmanI [200l . 
Preheating, the assumption of an elevated initial 
gas entropy at high red shift, was i ntroduced by 
lEvrard fc Henrvl ()1991l ) and lKaiseil (|1991f ) as a means to 
resolve the discrepancy in shape between the observed 
X-ray luminosity function and that expected from self- 
similar scaling of the cosmic mass function. The PH sim- 
ulation is tuned to match X-ray observations of clusters 
at redshift zero, particularly the lum inosity-temperature 
(L-T) relation (|Hartlev et al.ll2008f ). The preheating is 
achieved by boosting the entropy of every gas particle to 
200 keV cm^ at redshift z = 4. Although the preheat- 
ing dominates in the PH simulatio n, there is also cooling 
based on the coohng function of [Sutherland fc Dopital 
(|1993f l. Fewer than 2% of the baryons are converted to 
stars, however, and star formation is essentially halted 
by the preheating at z = 4. 

This simple model certainly does not capture all of the 
complex astrophysical effects associated with star and su- 
permassive black hole formation in clusters. The central 
entropy structure in local clusters is distinctly bimodal, 
with roughly half the population centered near the 200 
keV cm^ value used in the PH and the remainder cen- 
tered on an entropy level a factor of ten smaller. How- 
ever, the latter re flects potentially cyclical AGN feed- 
ing and feedback (IVoit et al.l 120021 : ISijacki et all l2007l : 
iPuchwein et al.l |2008[ ) that stronglv influences only a 
small mass fraction of the total cluster gas. Obser- 
vational evidence of u biquitous galactic winds at high 
and m oderate redshifts (jPettini et al.ll200ll : IWeiner et al.l 
|2008[) suggest that most of the heating of the ICM occurs 
at high redshift. Further support for fast feedback comes 
from red sequence galaxies extending t o redshift z = 1.4 
in th e Spitzer/IRAC Shallow Survey (|Eisenhardt et all 
|2008|) , the colors of which are consistent with passive evo- 
lution of a burst of star formation at redshift z ~ 3 — 4. 

2.2. Halo Catalog 

We identify halos in the simulation as spherical regions, 
centered on filtered density peaks, that encompass an 
average density of Ac/3c(-z), where Pc(,z) is the critical 
density of the universe. Both dark matter and baryons 
are included in the density measurement. We use ha- 
los identified at an overdensity of Ac = 200 for most of 
our analysis. Halo centers were identified with an A^*'* 
nearest neighbor approach, which approximates the local 
density by calculating the distance to the 32""^ nearest 
dark matter particle. The groupfinder begins with the 
dark matter particle with the highest local density, and 
works outward in radius particle by particle (including 
gas and dark matter) until the interior mean density is 
2QQpc{z). The algorithm then identifies the densest dark 
matter particle not already in a halo, and continues it- 
eratively until all overdense regions with more than 100 
particles have been identified. Overlapping halos are per- 
mitted; however, the center of mass of a halo may not be 
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Fig. 1. — Differential counts of halos versus total mass at redshift 
zero for the PH simulations (filled, black points), the GO simulation 
(open, red points), and the prediction from the Tinker mass func- 
tion (Tinker et al.„200S^ (solid, black line). 

in another halo. 

At redshift zero, in the PH simulation we have approxi- 
mately 220,000 halos with at least 100 particles, and 4474 
over a mass cut of M > 5 x lO^'^ h~^MQ. These num- 
bers are higher in the GO simulation, with approximately 
370,000 halos with at least 100 particles, and 5612 over 
the mass cut of M > 5x lO^^ft.~^Af0. In Figure[Tl we plot 
differential halo counts as a function of total mass from 
the GO and PH simulations, as well a s the predictio n from 
the Tinker mas s function (TMF) ^Tinker et all [20081. 
As discussed in IStanek et al.l (|2009D . preheating causes 
a decrease in total halo mass of up to ~ 15% relative 
to the GO treatment, with the largest effects at lower 
masses and higher redshifts. While the GO halo space 
density matches the TMF expectations well, the number 
of halos in the PH case is lower, especially at lower mass. 

Among the halos over the mass cut of M > 5 x 
10^'^ h~^MQ, we identify overlapping halos. For a pair 
of overlapping halos, we denote the less massive halo as 
a "satellite" and the more massive halo as a "primary" 
halo. For the rest of the analysis in this paper, we ex- 
clude the satellite halos. At redshift zero we are left with 
a sample of 4404 halos in the PH simulation and 5498 ha- 
los in the GO simulation. 

We repeat the halo finding exercise at all redshifts 
available for each model. In the case of PH, we employ 
a total of 63 outputs extending to a redshift of two. For 
the GO, we analyze only a subset of outputs at redshifts, 
z — 0, 0.5, 1.0 and 2.0. Unless otherwise noted, all of our 
analysis is for primary halos over the total mass limit of 
5 X 10^3 h'^ Mq. 

2.3. Bulk Halo Properties 

With the primary halo samples identified, we calcu- 
late bulk properties for them that we roughly classify 
into "structural" and "observable" categories. The for- 
mer includes dark matter velocity dispersion, ICM mass 
fraction, gas mass-weighted temperature, and halo con- 
centration while the latter includes X-ray luminosity 
and spectroscopic-like temperature, thermal Sunyaev- 
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Zel'dovich effect, and a dimensionless ICM emission mea- 
sure. 

As the MGS simulations are SPH treatments of the 
gas, integrals over volume map to summations over all 
particles, via fdVp"^ T,imipf~^ . We consider two 
measures of ICM temperature. First, we consider the 
mass-weighted temperature. As GADGET-2 is a La- 
grangian simulation with equal mass gas particles, the 
mass-weighted temperature is simply the average tem- 
perature of the particles in the halo: 



T„ 



1 

M 



dVpT 



1 ^ 



(1) 



We also calculate t he spectros c opic-l ike temperature, 
Tsi, as defined in Mazzo tta et~aLl ()2004[ ). 



1-1/2 



-3/2 ' 



(2) 



with a = 0.75. The spectroscopic-like temperature offers 
a good match to the temperature derived from a one- 
component fit to an X-ray spectrum, but is far simpler 
to compute. 
We calculate X-ray luminosities. 



dVp^K[T) 



N 

E 



AA(r,), 



(3) 



using MEKAL tables assuming fixed 0.3 solar metallicity 
to calculate A(T). We do this in energy bands of 0.7-2.0 
keV, 0.7-5.0 keV, and 0.7-7.0 keV, and also compute a 
bolometric luminosity, Lboi using a wide photon energy 
range of [0.1-40.0] keV. These tables include both con- 
tinuum and emission lines derived from an a ssumption 
of coUisional ionization equilibrium ()Mewe et al. 2003). 

We also calculate the global thermal Sunyacv- 
Zel'dovich si gnal, Y , parameter, f o Uowing the conventio n 
presented inlda Silva et al.l (l2000D:r Springel et all (|2001[): 
iKravtsov et al.l (|2006D . 



Y = 



dVupTp 



gas 



N 

meC^ A rrip J 

' (4) 

With fcs Boltzmann's constant, ot the Thomson cross- 
section, me the electron mass, and rrip the proton mass, 
and A the halo comoving area. The above expression 
yields Y in units of h~'^ Mpc^. 

2.4. Radial Profile Measures 

In addition to bulk cluster properties, we include two 
measures of radial structure in our analysis: halo concen- 
tration of the total mass and a dimensionless emission 
measure for the hot gas. 

We measure the ha lo concentration, c, after fitting an 
NFW density profile (|Navarro et al.lll997D . 

Pcrit {r/rs){l + r/rs)^' 

to the radially-binned, total mass density profile of each 
halo. Following NFW, we take the concentration, c, to 
be defined at a density threshold of 200 times the criti- 
cal density, c = r2oo/rs. Halo concentration is a good 



proxy for halo fo rmation epoch (jWechsler et al.l 120021 : 
iBusha et al.ll2007l ). We show below that baryon loss in 
the PH case has a non-neglgible effect on halo concentra- 
tions. 

We also measure a dimensionless emission measure, or 
clumping factor, Q, that measures the contribution of gas 
density structure to the halo luminosity. Using a scaled 
radius, y = r/r2oo, we write the radial gas density profile 
of a given halo as 

p{yr20o) = /icM(200pc)g(2/), (6) 

where /icm = Mgas{< f'20o)A^20o is the halo's ICM mass 
fraction within r2oo- With the overall ICM mass frac- 
tion factored out, g{y) becomes a dimensionless structure 

function normalized by 3 dyy^g{y) = 1. 

The second moment of g(y) defines the dimensionless 
emission measure 



Q(T) = (3/47r) / dhjg\y). 



(7) 



This definition allows the X-ray luminosity scaling to be 
written as 

2 A 



L cx pcM2QoA{T) /icm' 



(8) 



We show below how lower values of both /icm and Q for 
the PH case drive X-ray luminosities down by an order of 
magnitude at lO^'' Mq relative to the GO treatment. 

3. MEAN SCALING RELATIONS 

To characterize the mean behavior of halo properties, 
we perform a linear least squares fit to the natural log of 
the i"^ signal as a function of mass and redshift, using 
the form 

{si){p,a) Sisiia) + ai{a)p, (9) 

where = InS^ and p = ln(M/ H-^IO^'^Mq). Brackets 
represent averaging in narrow mass bins at fixed epoch, 
so ai{a) is the slope and 3^^14(0) the normalization at 
lOi-* h-^MQ of the i^^ signal at redshift, z = a'^ - 1. In 
the PH case, one expects the entropy floor to introduc e 
curvature into the scaling relations (|Voit et al.l I2002D . 
and so we extend this model to a quadratic with respect 
to p for the ICM mass fraction and related measures for 
this case. 

At each redshift, the mass scalings are derived from 
the sample of primary halos with M > 5 x 10^'^ M©. 
We then fit the evolution of the normalization, presented 
as a shift relative to the present value, 

Asis4.ia) = Sisiia) - Sj,i4(l), (10) 

to the form 

As,,i4(a) = /31n(S(a)), (11) 

where E{a) = H{a)/Ho. For the PH model, this form is 
a poor fit to /icM-related quantities. We find a better fit 
in terms of a quadratic in In a, 

Asi,i4(a) = eo + ei ln(a) + e2{\n{a)f, (12) 

In this section, we present mean scaling relatio ns and 
discu ss deviations from self-similar expectations f Kaiseil 
Il986l ) for structural ( §3.ip . and observable ( §3.2|1 proper- 
ties. We then demonstrate good agreement between X- 
ray observations and the low-redshift PH model results. 
Deviations from mean scaling behavior, in the form of 
a signal covariance matrix, ((si — {si)){sj — (sj))), are 
presented in 21 
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Fig. 2. — Scaling relations at redshift zero (left) and redshift evolution of fit parameters (right) for ctdMi Tm, /iCM> and the NFW 
concentration c (top to bottom). The PH model data are black, filled points and the GO are red, open points. Subsamples of halos are shown 
in the left panels, drawn in a manner that shows ~ 10% of the overall population with nearly uniform mass sampling. Solid lines show the 
best-fit sc alin gs derived from the entire sample. The right-hand pan els p resent the evolution of the slope and the shifts in normalization, 
equation HIOII . Solid lines show fits to log-linear behavior, equation mil l, except ficm which is quadratic in InM, while the dotted lines 
give the self-similar predictions listed in Table \3\ 
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3.1. Structural Quantities 

The left panels of Figure [5] present scaling relations as 
a function of mass at z = for sub-samples of the GO 
and PH halo samples. Four structural measures are pre- 
sented: dark matter velocity dispersion, ctdm; intraclus- 
ter gas mass fraction, /icm; mass- weighted gas tempera- 
ture, Tm', and NFW concentration, c. Best- fit parameters 
to the mass scaling, equation ([9]), are presented in Table 
[T] The right panels show redshift evolution of the slopes 
and the shifts in normalization, equation (jlOp . for each 
signal. Error bars in the fit parameters are derived from 
bootstrapping resampling of the samples. In general, the 
uncertainties on the best fit parameters are very small, 
~ 0.1% at redshift z = 0, and are much smaller than the 
intrinsic scatter at fixed mass about the median power 
law relation. The errors grow larger at higher redshifts; 
the mass-limited sample size drops below 100 at z — 1.8 
for the PH case. Fits to the redshift evolution of the nor- 
malization, equation pip are presented in Table [3] 

Dark matter velocity dispersion. The velocity disper- 
sion of the dark matter particles is a fundamental mea- 
sure of the virial state of a halo. Observationally, the 
galaxy velocity dispersion tracks ICM temperature in a 
manner consistent with virial expectations, but the pos- 
sibility of a ^ 10% bias relative the dark matter is still 
allowed (|Biviano fc Katgertl[200llBecker etalMOTh . In 
both simulations. Figure [2^ shows that the scaling with 
mass is slightly shallower than the self-similar predic- 
tion, (TDM (X M^/^. The dark matter velocity disper- 
sion at 10^^ Mq is in good agreement with the 
1082.9 ± 4.0 kms~^ value derived from a suite of N- 
b ody simulations bv lEvrard et al.l ([2008). As discussed 
in lEvrard et al.l (l2008i r a slight suppression of the slope 
is expected from finite particle resolution when the low 
mass halo cutoff corresponds to a few thousand particles, 
as is the case here. The fact that the PH and GO slopes 
are similar (0.341 and 0.345) indicates that the difference 
in gas physics treatments is not responsible for the shift. 
In both PH and GO cases, the slope remains constant to 
high redshift, as seen in the upper right panel of Figure 
[21 The evolution of the normalization, (Tdm,14 shown 
in Figure [2j3, deviates slightly from the self-similar pre- 
diction, ctdm oc [E{z)M]^/^ , with the velocity dispersion 
being 1 — 2% higher at z = 1 (values of (3 are given in 
Table [3]) . Despite strong baryon content differences dis- 
cussed below, the dark matter virial scaling under the PH 
and GO treatments remain remarkably consistent. 

Mass-weighted temperature. The mass-weighed tem- 
perature, Tm, a useful probe of the hydrodynamic 
state of the ICM, is known to have small scatter 
with respect to rnass in simulations (lEvrard et al.lll996t 
iBrvan fc NormanI Il998t iBorgani et al.ll2004( l. Figure [2fc 
shows that the slope of the GO scaling relation agrees 
with the self-similar scaling, T ~ M^/^, at the few per- 
cent level, whereas the PH slope is significantly less steep, 
a = 0.559 ± 0.002. As all gas particles receive the same 
entropy boost at z ~ 4 in PH , the resultant fractional 
increase in the characteristic entropy of a halo is larger 
for low-mass halos than for high-mass halos. The ele- 
vated initial entropy drives the tilt as well as a ~ 35% 
higher normalization of the Tm — M relation in the PH 
simulation. But the effect diminishes at higher masses; 
the PH halos at 10^^ h'^ Mq are only 10% hotter than 
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Fig. 3. — The normalized bulk kinetic energy in gas turbulent 
motions, /3gas = o"gas^/(fcTm/Atnip), is shown as a function of mass 
for sub-samples of the PH (filled, black points) and GO (open, red 
points) halos. Lines show best-fit mean scalings. 



their GO counterparts. 

The slopes of the Tm — M relation do not evolve 
strongly with redshift in either model. The normaliza- 
tions, however, do not follow the self-similar expecta- 
tion, with the GO simulation lying ~ 10% low and the PH 
case 5% high at z — 1 relative to Tm ^ E{a)'^^^ scaling 
(|Kaiseilll98 6^. The self-similar prediction is for halos in 
perfect hydrostatic equilibrium, and mergers are known 
to drive deviations from hydrostatic equilibrium. How- 
ever, the total kinetic energy of the gas — the sum of 
thermal plus bulk kinetic, or turbulent components — 
should scale more closely to self-similar expectations. 

We define a dimensionless measure of the kinetic en- 
ergy content of the gas, /?gas = c^gas'^ /{kTm/ fJ-rUp), where 
/i is the mean molecular weight of the gas. Figure [31 
shows values of /3gas for a subsample of halos at z = 0. 
The values differ strongly in the two simulations, with 
mean values at lO" h'^ Mq of 0.16 in the GO simula- 
tion and 0.06 in the PH model. The positive scaling with 
mass reflects the later formation epoch of high-mass ha- 
los. The next generation of X-ray telescopes, such as 
IXO^, will have high-resolution spectroscopy capable of 
measuring figas, and thus be able to discriminate between 
the model predictions in Figure [3l 

We measure the total gas kinetic energy, 



Etnt — 



kTm 
limp 



'-''gas J 



(13) 



and fit this to the standard form, equation ([9]). Figure 
m shows the redshift evolution of the total gas kinetic 
energy in both models. The GO simulation respects the 
self-similar scaling, E^tot ^ [ME{a)]'^^^ at the few percent 
level to redshifts z = 1. At a given mass, the velocity dis- 
persion contribution to the total kinetic energy increases 
with redshift, and this enhanced turbulence is responsi- 
ble for driving the mass- weighted temperature away from 
self-similar scaling, as seen in Figure [2l For the PH case, 

http://ixo.gsfc.nasa.gov/ 
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TABLE 1 

Redshift Zero Mass Scalings ^ 



PH Simulation 



GO Simulation 



Signal 


Si4 


a 


Si4 


a 




6.1990 ±0.0006 


0.341 ±0.001 


6.2037 ± 0.0006 


0.345 ±0.001 


kT ^ 


0.5121 ±0.0008 


0.559 ± 0.002 


0.209 ± 0.001 


0.650 ±0.002 




0.605 ±0.001 


0.576 ± 0.002 


0.163 ±0.003 


0.576 ±0.005 


yb 


-12.810 ±0.002 


1.825 ±0.003 


-12.642 ±0.002 


1.651 ±0.003 




-1.653 ±0.003 


1.868 ±0.006 


0.622 ± 0.005 


1.079 ± 0.008 


NFW c 


0.966 ± 0.005 


0.071 ±0.007 


1.181 ±0.004 


-0.002 ± 0.006 


Q 


0.461 ± 0.002 


0.161 ±0.003 


1.219 ±0.001 


0.030 ± 0.003 



Fit Parameters to equation JgJ with uncertainties from Monte Carlo re- 
sampling. 

Units: (TDM (kms-i); fcTm or kT^i (keV); Lboi (10'''' erg s^^ cm^^). 
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TABLE 2 

Redshift Zero PH Quadratic Mass Scalings 



Fig. 4. — The evolution of the normalization and slope of the 
total kinetic energy of the gas in the PH (filled, black points) and 
the GO (open, red points) simulations. The dotted line is the self- 
similar prediction and the solid line is the measured evolution in 
the PH simulation. 

the degree of turbulence is much lower, and the total en- 
ergy remains higher at high redshift due to the influence 
of the initial entropy injection. 

Baryon fraction. Figure [2fc,f shows the scaling of the 
hot gas fraction, /icm, with mass at z = along with the 
redshift evolution of the slope and normalization. For 
reasons discussed below, we show evolution in the PH 
normalization and slope at a mass of 5 x 10"'^'* M© as 
well as the fiducial mass of 10^'' M©. 

At redshift zero, the distribution of /icm values within 
^200 in the GO simulation is simple: halos have a slightly 
depleted baryon fraction, with mean /icm — 0.90 , and 
a dispersion of 0.04. There is no trend with mass, as 
the slope is 0.001 ± 0.019, a result consistent with pre- 
vious SPH simulati ons done at simila r re solution and 
evalua ted at r2oo bv lCrain et al.l (|2007f) and lEttori et al.l 
(I2006D. The simulation of the MareNostrum universe 
iGottlober fc Yeped ()2007t ) has twice our mass resolution, 
and a baryon fraction of /icm = 0.92 at a virial radius, 
Tmri that encompasses a mean density of ^ 100 times 
the critical density. This small but significant increase is 
consistent with the radial trend seen in our simulation. 
O ur value of /icm i s lowe r than the 0.97 value measured 
bv lKravtsov et al] (|2005f) at in their adiabatic AMR 



Signal 


Sl4 




«2 


/icm 


-2.270 ±0.002 


0.310 ±0.009 


-0.0661 ±0.0061 


yb 


-12.790 ±0.005 


1.864 ±0.018 


-0.0586 ±0.0017 


T b 


-1.641 ±0.012 


1.892 ±0.043 


-0.0351 ±0.0042 



^ Fits of the PH signal-mass relations at 2 = to a quadratic, 
ln5 = si4 + ailnAf + a2(lnM)^, with mass in units of 
IQl'' h-^ Mo 



10. 



Units: Y (h"^ Mpc^); Ltoi (10"" erg s" 



TABLE 3 
PH Normalization 
Evolution in ln(i?(a)) 



Signal 


/3 


Self-Similar 


o"DM 


0.34 


1/3 




0.76 


2/3 


Ta 


0.73 


2/3 


Lhoi 


1.39 


7/3 


/icm 


-0.44 





Y 


0.33 


2/3 


Q 


-0.13 


N/A 


c 


-0.68 


N/A 



^ Equation [TT] 

(Adaptive Mesh Refinement) code. That study included 
a comparison of SPH and AMR simulations evolved un- 
der gravity only, and they note a statistically-significant 
offset of ~ 5%, with the AMR simulation having higher 
baryon fraction. 

As has been shown in previous simulations 
(Bial elt et al. 2001; Borg ani et aI]|2001l:lMuanwong et all 
120061 : lYounger fc BrvanI |2007[ ). preheating has a dra- 
matic effect on the hot gas fraction of massive halos. 
Low-mass halos can lose half of their baryons within 
^200 while the most massive halos are depleted by only 
~ 10% relative to the GO case. A power-law form is a 
poor representation of the PH mean behavior of /icm 
with mass, so we extend the model to a quadratic in 
In M and give best-fit parameters in Table [51 We present 
a comparison with observed gas fractions in i i3. 31 below. 

In the GO simulation, the slope of the gas fraction 
remains consistent with zero at all redshifts, while the 
mean gas fraction increases slightly, rising from /icm = 
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TABLE 4 
PH Normalization Evolution in ln(a) 



Signal 








^2 


/iCM 


-7.64 X 10" 


-4 


0.0143 


-0.371 


Y 


-1.49 X 10" 


-3 


-0.308 


-0.0141 




-9.44 X 10" 


-4 


-0.786 


0.477 



^ Equation [T2l 



0.90 at z = to 0.93 at z = 1. The latter effect is 
consistent with the findings of lGottlober fc YepesI ()2007l ) 
in the Marenostrum simulation, where the mean baryon 
fraction increases from 0.92 at z = to /icm — 0.94 at 
z = 1. While the mechanism responsible for this slight 
drift is not fully understood, energy transfer from the 
dark matter to the gas during mergers m ay play a role 
(|Pearce et al.|[l99l iMcCarthv et al.|[2007l) . 

The redshift evolution of the baryon fraction in the PH 
simulation depends on mass scale. Gas expelled from 
low-mass halos by preheating at z ~ 4 can accrete back 
onto descendant halos at later stages in the merger hi- 
erarchy, in a manner that depends on the characteristic 
entropy of the later halo. At the fiducial 10^"* Mq 
normalization scale, halos at z = 1 have a ~ 20% lower 
gas fraction compared to z — 0, and the local slope of the 
/icm — M relation is also steeper at higher redshift. The 
effects of preheating dominate over the universal expan- 
sion at the 10^'^h~^MQ mass scale; hence the evolution of 
the normalization cannot be simply described as a power 
of E{z). For the baryon fraction, and for other halo pa- 
rameters which depend strongly on the baryon fraction, 
we fit the evolution as a quadratic in ln(a), with the 
best-fit parameters presented in Table ID 

At a higher mass scale of 5x lO^^/i~^M0, the effects are 
milder, and the local slope is very close to the slope of the 
GO model. The most massive halos are less affected by 
preheating, and their baryon fractions are ther efore more 
appropriate to use as indicators of cosmology (|Pen et al.l 
120031 (Alien et al.ll2004L l2008l ). We note that the mean 
gas fraction shift at 5 x 10^'' Mq in the PH model 
is comparable to t he 20 percen t uniform prior on gas 
fraction applied by lAUen et al.l (.2008) in their analysis 
of a Chandra sample of 42 clusters with kT > 5 keV 
extending to z = 1.1. The cosmological constraints from 
that work would thus not be strongly affected if a PH 
model prior on /icm behavior were imposed. 

NFW concentration. Figure [2|;,h show the behavior of 
the NFW concentration, c, derived from fits to the to- 
tal mass density profile of each halo. At 10^^ M©, 
the mean concentration in the GO simulation, c = 3.3, is 
substantially lower than the value of 5.2 found for the 
original Millenniu m Simulat ion w ith only dark matter 
(|Neto et al.l 120071 : iGao et al.if2008l) . The mean concen- 
tration in the PH simulation, c = 2.6, is lower than the 
GO value. As discussed above, the PH halos have a lower 
baryon fraction than GO halos, and the increasing baryon 
loss at low masses has the effect of tilting the c — M re- 
lation so that the slope is positive rather than weakly 
negative. 

3.2. X-ray and SZ Signals 

We now turn to common bulk observed properties of 
clusters: the SZ decrement, Y , and X-ray spectroscopic- 



like temperature, Tgi, luminosity, iboi, and the emission 
measure Q. 

Sunyaev-Zeldovich decrement. From Figure [2t,e, we 
know that preheated halos at fixed mass have, on aver- 
age, lower ICM mass fractions and higher mass-weighted 
temperatures than their GO counterparts. Since the in- 
tegrated thermal SZ decrement, Y, is a product of these 
two measures, we can anticipate some degree of cance- 
lation coming from this opposing behavior. Figure [5^ 
shows that this cancellation is quite close to exact above 
a mass scale of 2 x 10^** Mq. 

The slope in the GO model is very close to the self- 
similar expectation of 5/3, but curvature in /icm tilts 
the PH relation from a local slope of —1.8 to —1.6 across 
the mass range IO^'^-IO^^H-^Mq. As with the ficM-M 
relation, we fit 1^ to a quadratic in In M for the PH case 
in order to account for the curvature and to get the best 
possible measure of the intrinsic scatter about the mean 
relation. Parameters are given in Table [2l 

As discussed further in 21 the dispersion about the 
intrinsic, mean Y — M relation is 12 ± 1 percent in both 
models. The tight mass scaling is independent of clus - 
ter dynamical state, as noted bv iKravtsov et al.] (|2006l ). 
making the SZ thermal signal ideal for cluster detec- 
tion. It must be remembered that this scatter is en- 
larged by line-of-sight projections, and t hat the magni- 
tude o f this effect w i ll be model-dependentlSpringel et al.l 
((200lh : lShaw et all (|200l ): iHallman et all (|2007D . 

While the Y — M normalization at high masses is not 
sensitive to cluster physics in our models, previously 
ART simulations with cooling, star formation and su- 
pernova feedback (CSF) displ ay a normaliz ation drop by 
25% relative to the GO case iNagail (|2006l ). At a basic 
level, the different behavior between the two studies sim- 
ply refiects the fact that the thermal pressure support, 
— VP/p, is insensitive to a multiplicative shift in gas den- 
sity normalization. Our models force all baryons into the 
hot phase, while a CSF treatm ent allows som e fraction of 
baryons — 40% in the case of lNagail ()2006f ) — to reside 
in galactic sinks of stars and cold gas. Since observations 
indicate that the ratio of stellar to hot gas mass declines 
with increasing mass, fr om values near 0.5 at 10^^h~^MQ 
to 0.1 at 10^5 h-^ Mq (|Giodini et al.ll2009( ). one would 
anticipate that the observed Y — M relation to be steeper 
than self-similar. A recent X-ray based estimate of the 
Y ~M relation from a Chandra archival sample of groups 
and clusters finds a slope of 1.75 at rsnn o ver the mass 
range lO^^ - lO^^ /j-i M© (|Sun et al.ll20Ml . 

Figure compares the redshift evolution of the Y — 
M relation in the two simulations. The slope in the PH 
model steepens at higher redshift, and the normalization 
drifts below self-similar expectations by ~ 10% at z = 
1. The evolution at high mass (~ 5 x 10"^* Mq) 
in the PH model is similar to that of the GO model, in 
both slope and normalization, since redshift z ~ 1. At 
our chosen normalization mass of ~ 10^* Mq, the 
evolution departs from self-similar at z > 0.5, driven 
by the evolution of the baryon fraction at that scale. 
Like the baryon fraction, we fit the evolution of y to a 
quadratic in ln(a), with the best fit presented in Table 
m We note that the PH behavior is mildly in conflict 
with the CSF model evolution of Nagail (|2006f ). which 
displayed consistency with self-similar evolution at the 
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~ 20% level. With only 11 halos, that study could not 
address statistical differences at the level we do here. 

Spectroscopic-like temperature. We also consider the 
spectroscopic-lik e temperature, Tgi, an analytic prescrip- 
tion derived by 'Mazz otta et al.l (|2004l ) to approximate 
observed X-ray spectral temperatures. We find that Tgi 
agrees well with the X-ray temperatures d erived from 
spect ral fits to X-MAS2 mock observations ll Rasia et al.l 
|2005[) . As seen in Figure [SJ;,d, the slopes in the PH and GO 
treatments agree well, a = 0.57, but at fixed mass the PH 
halos are ~ 40% hotter than the GO halos. The distribu- 
tion of the gas in the GO mod els contains cool, low entrop y 
cores of accreted sub-halos (jMathiesen fc Evrardll200ll ). 
and these cool, dense clumps pull down the Tgi measure 
relative to the PH simulation. The cores in the latter case 
have been effectively erased by the preheating. 

The redshift evolution of Tsj. 14(a) in the PH simulation 
is very close to self-similar (cx [i5(a)]^/^), and tracks well 
the mass- weighted behavior shown above (see Table [3]). 
The evolution in the GO simulation departs dramatically 
from self-similar behavior at low redshifts. 

Bolometric luminosity and Emission measure. The dif- 
ferent behaviors of the baryon fraction and temperature 
in the PH and GO simulations drive X-ray luminosity dif- 
ferences, but the gas clumping, or emission measure, also 
plays a significant role. Figure [5]3,g show the z = mass 
scalings of the bolometric luminosity, L^oi, and dimen- 
sionless emission measure, Q. Because of the influence 
of line emission, the GO scaling is less steep than the self- 
similar slope of L ^ M**/^, and the emission measure is 
nearly constant with mass at a value ~ 3.5. In the PH 
case, the clumping factor is smaller by a factor of two 
and displays a significant trend with mass. Combined 
with the /icM behavior, the result is a suppression of 
Lboi in the PH case by a factor of 10 at 10^'' Mq, and 
a steepening of the slope in mass to 1.87, from 1.08 in 
the GO case. 

Although the PH halos have a higher temperature at 
fixed mass, this effect is dwarfed by the decreases in /icm 
and Q between the PH and GO treatments. The lower nor- 
malization of the Q — M relation reflects a shallower gas 
density profile in the PH case. In turn, the lower central 
gas densities contribute to lowering the overall mass pro- 
file, driving the shift in concentrations to lower values for 
the PH halos discussed above. As in the Y — M relation, 
the curvature in the /icm — M relation in the PH model 
drives curvature in the Lboi — M relation. iHartlev et al.l 
([2OO8I) see this curvature in the L — T relation in the 
PH model, as well as in a large, local observed sample of 
galaxy clusters. We fit the Lboi — M relation in the PH 
model to a quadratic in InM, presenting the best fit in 
Tabled 

The evolution of the normalization in PH is not a per- 
fect power of E{a), due to the complicated evolution of 
/icm with redshift. We fit the evolution of ^14(0) to a 
quadratic in ln(a) in the PH simulation, and note that it 
is weaker than the GO evolution. However, as shown in 
Figure [5]:, this evolution is a function of mass. The larger 
solid points, for halos at 10^^h~^MQ, show a weaker evo- 
lution than the evolution of the halos at 5 x 10^** Mq , 
shown by the smaller solid points. The latter halos evolve 
similarly to the GO halos, illustrating that the most mas- 
sive halos in the PH model are similar in structure and his- 



tory to the GO population. The evolution of the iboi ~ M 
relation in the PH simulation is driven mainly by the red- 
shift evolution of the hot gas fraction, with Q contribut- 
ing 10% of the decrease at z = 1. 

3.3. Comparison to Observations 

As shown bv IHartlev et al.l ()2008l ). the PH simulation 
offers a good match to the core-excised L — T relation of 
local clusters. Here, we briefly explore the level of agree- 
ment between the models and observations in this and 
other scaling relations. While not wishing to oversell the 
simple physical treatment of the preheated model, which 
is undoubtedly wrong in detail, we show below that it 
reproduces several scalings with quite high fidelity. Pre- 
heating appears to be a useful effective model. It is 
important to remember that precise comparison of ob- 
servations and simulation expectations requires careful 
modeling of survey selection and projection effects, and 
we do not treat these effects here. 

Figure [6] compares the Lboi — Tgi relations with the 
core-excised iboi — Tx measur ements of the loca l, rep- 
resentative REXCESS survey (jPratt et al.ll2008D . The 
slope of the observed relation is somewhat shallower, and 
its scatter somewhat larger, than the PH model predic- 
tions, but these differences are at the level of a few tens 
of percent in luminosity, or less than ten percent in tem- 
perature. The small scatter in the core-excised Lboi — T 
relation for REXCESS clusters (and PH halos) indicates 
that local galaxy clusters are we ll-behaved outside of the 
core ("Neu mann fc Arnaudll2001h . 

At higher redshift, we consider the CCCP clusters, a 
subset of the 400 square d egree survey which h as been 
followed- up with Chandra (jVikhlinin et al.ll2008f) . These 
clusters range from approximately 0.3 < z < 0.8, so at 
fixed mass we scale, in a self-similar manner, the ob- 
served luminosities and temperatures to z = 0.5 for com- 
parison with the model Lboi — ?si relations. The latter 
are measure d within = 500 to b e consistent with the 
treatment of IVikhlinin et al.l ()2008D . The comparison is 
presented In Figure [T] The agreement is good, but the 
measurement errors are larger than the REXCESS sam- 
ple. The scatter in the CCCP sample is larger, and this 
may be due in part to the fact that cores have not been 
excised in the luminosity measurements of this sample. 

Figure [8] compares the mass scaling of ICM mass frac- 
tions at redshift zero in the models to XMM measure- 
ment s for local clusters ()Arnaud et al.l l2007l : ISun et all 
l2009l) . Values are measured within Ac = 500. To esti- 
mate observed cluster masses, hydrostatic estimates that 
include a radial temperature gradient are employed. We 
note that Arnaud data agr ee with other observational de- 
term inations at hi gh mass (lVikhlininll2(MlGiodini et al.l 
I2009D . while the iSun et all (|2009f) data extend to lower 
mass, ~ lO^'^ Mq. The PH model matches the ob- 
served /icm — M relation well, with hot gas fractions a 
factor two less than the cosmic ratio, ftb/^m at 10^^ and 
a strongly increasing trend toward higher masses. Larger 
observed samples are needed to test the curvature and 
degree of scatter. 

In Figure ^ we compare the T s \ — M relation of the 
models to that from lArnaud et al.l ()2007f ). The GO model, 
with its cool sub-halo cores, is strongly offset from the 
data. The PH model relation is much closer, with a sim- 
ilar slope and scatter. There is, however, a consistent 
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Fig. 5. — Scaling relations at redshift zero (left) and redshift evolution of fit parameters (right) for Y, T^i, Lboli a^nd Q (top to bottom). 
Point and line styles are identical to Figure |2] 
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Fig. 6. — The redshift zero Lbol — Tsi relation for GO (open, 
red points), PH (filled, black points), measured within Ac = 500, 
and the core-excised I/bol ~ T relation from the REXCESS survey 
(large, blue points) llPratt et al.ll2008r) . 




1 10 
kT,i [keV] 

Fig. 7. — The redshift z = 0.5 ip. 5-2.0] - ^sl relation for PH 
(filled, black points), measur ed within Ac = 5 00, an d the CCCP 
L[Q 2—2.0] ~ relation from IVikhlinin et al.l 1 12008) (large, blue 
points), scaled to redshift z = 0.5. 

offset of ~ 15% in mass toward lower values in the ob- 
served sample. The magnitude of this offset is consistent 
with the level of expected bias from hydrostatic mass es- 
timates, which simulations show tend to underestimate 
true masses by a pproximately 20% (jRasia et al.l 120061 : 
iNagai et al.ll2007t ). 

Overall, the bulk X-ray properties of the PH simula- 
tion match the local scaling relations quite well. This is 
particularly true after excising the core from the obser- 
vations, and after considering the mass bias introduced 
by hydrostatic mass estimates. Although the physics of 
the PH model is simple, it appears to provide useful rep- 
resentation for the behavior of the bulk of the hot ICM 
lying outside the cool core regions. 
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Fig. 8. — The redshift zero /icM — M relation for GO (open, red 
points), and PH (filled, black points), measured within Ac = 500. 
The observations plotted are the data f rom Arnaud et al. (200^ 
(large, blue points) and from lSun et al.l II200&) (green crosses). 
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Fig. 9. — The redshift zero T^i — M relation for GO (open, red 
points), PH (fi lled, black points), m easured within Ac = 500, and 
the data from lArnaud et al.l II2007I) (large, blue points). 

In this section, we explore the second moment of the 
halo scaling relations. Understanding the variance at 
fixed mass is necessary for calculating the mass selection 
properties of signal-limited samples, and survey counts 
typically constrain only a line ar combination of signal 
normalization and variance (St anek et aLll2006l ). In ad- 
dition, the signal covariance determines the precise struc- 
ture of scaling relations in signal-selected samples. A 
worked example of the L — T rel ation expec t ed for X-ray 
flux-limited samples is given bv lNord et aP (|2f)f)8h . 

After defining terms, we begin by presenting the covari- 
ance of signals at fixed mass at the present epoch, then 
demonstrate that redshift evolution in most elements is 
weak. We close with analysis of the mass selection prop- 
erties of signal pairs. 



4. COVARIANCE OF BULK PROPERTIES 



4.1. Signal Covariance Matrix 
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For a set of halos of mass at expansion epoch a, we 
define a symmetric covariance matrix, with elements 

*u = {{si - Si{fi,a)){sj ~ Sj{fi,a))) (14) 

where the mean values are determined by equation ^ 
and the brackets represent an ensemble average. The 
J*'' diagonal element of the covariance matrix is the j*'' 
signal variance, cr| = V^jj. We present these measures 
along with the correlation matrix, Cxy = '^xy/{<^x'^y), 
which expresses the covariance in normalized terms. 

Figure [10] presents a graphical representation of the 
data comprising the covariance matrix for a subset of 
signals at redshift zero. The diagonal panels show the 
distribution of signal deviations (in the natural log of 
the measured signal) about the mean for the PH (shaded) 
and GO (lines) treatments. Panels off the diagonal plot 
the normalized deviations {{sj — Sj)/aj vs. {si — Sij/ui) 
for signal pairs, with the lower and upper triangles show- 
ing PH and GO cases, respectively. The orientation and 
spread of halos in each panel determines the correlation 
coefficient. For instance, the tight ellipse formed by the 
population in the Y — /icm panel indicates a high corre- 
lation coefficient between this pair of signals. 

The 2 = values of the correlation coefficients, along 
with uncertainties from bootstrap resampling, are pre- 
sented in Table[6l As in the mean fit parameters, typical 
uncertainties are on the order of ~ 1%. Before explor- 
ing the off-diagonal terms, we first examine the variance 
and distribution function shape of individual signal de- 
viations. 

4.2. Signal Variance at Fixed Mass 

The assumption of log-normal variance is a common el- 
ement of the likelihood analysis used in cluster cosmology 
studies (see references in §1). Both the absolute variance 
and the full distribution shape affect the expected signal 
counts derived from mass function convolution. 

In Table [HI we list root mean square deviations for the 
full set of signals at z = 0. To test the log-normal expec- 
tation, we list the normalized deviates at which each sig- 
nal's ranked distribution reaches fixed percentile values, 
taken to be ±2ct, iblcr and median/mean of a Gaussian 
distribution. The difference between the listed values 
and their integer counterparts is a measure of the de- 
gree of local deformation from Gaussian in the frequency 
distribution. 

While deviations from Gaussianity are apparent in es- 
sentially all measures, the typical percentile shifts in the 
PH model are only a few percent. Exceptions are a sig- 
nificant positive skew in the dark matter velocity dis- 
persion, with median location of — 0.12cr and shifts in 
the ±2(7 Gaussian tails to — l.Tcr and 2.5(7, respectively. 
Most distributions are slightly leptokurtic, especially the 
mass weighted temperature. Worth noting is the fact 
that the shapes of two important cluster selection ob- 
servables, Y and L\^oh do not deviate by more than 0.1 
from the Gaussian expectations. 

Under the GO treatment, the shape of the dark matter 
velocity dispersion is the same as in the PH case, showing 
positive skew at the 10% level. However, the shapes of 
the hot gas properties of GO halos generally differ, to a 
slight degree, from the PH shapes; all measures tend to 
be slightly skew negative. The differences can be sub- 



tle, as close inspection of the Y histograms in Figure [TOl 
confirms. 

Comparing the GO and PH distributions, we conclude 
that a log-normal approximation is a fairly accurate de- 
scription, but calculations demanding better than ~ 10% 
precision in shape will require an expanded treatment, 
either via direct Monte Carlo from simulations, or from 
analytic extension s, such as an Edgeworth series in Her- 
mite polynomials (jKofman et al.ll993tlShaw et al.l2 0091. 
We leave detailed analysis of the physical mechanisms 
driving these distribution shapes to future investigations. 
In the case of ctdm, we have preliminary evidence that 
mergers drive the positive skew tail of the deviations; 
nearly all halos with a deviation of > Scr have undergone 
a merger since redshift z = 0.2. 

The amplitudes of the scatter in the set of signals range 
from a low of 0.036 for /icm to a high of 0.28 for Lboi, 
both under the GO treatment. For the PH case, iboi 
is highest, at 0.19 while ctdm is lowest at 0.042. The 
latter value matches the GO case, and both agree with 
the scatter derived fro m the ensemble value presented 
bv lEvrard et all (|2008[ ). While the dark matter virial 
scaling is robust to simp le physical treatments for the 
baryons, ()Lau et al.ll2009D find that strongly dissipative 
baryon physics depresses the slope of the ctdm — M re- 
lation by introducing mass- and redshift-dependent in- 
creases in halo velocity dispersion. 

By raising the halo sound speed, which drives the shock 
radius to larger v alues and lowers the Mach number of 
infalling material (jVoit et al.ll2002l ). preheating leads to 
more thermally regular halo gas, with smaller scatter 
in mass-weighted temperature, T^, compared to the GO 
simulation. Naively, in a virialized cluster dominated by 
only gravitational effects, one expects (Tdm^ ^ 7m, and 
2fTcrpj^ — "'Tm • The values for the GO simulation are close 
to this, with aT„^ 2.5(T(7dm- However, the scatter in 
Trn is much lower in the PH simulation, only 1.4fTcrj2j^j , or 
5.8%. 

The scatter in /icm shows the reverse behavior. In the 
GO simulation, the scatter is very small, 3.6%, consistent 
with other gravity-only SPH simulations (Ettori ct al] 
l2006t ICrain et aLlboOTn . When only gravity drives gas 
thermodynamics, the gas distribution traces the dark 
matter distribution very well. The PH model scatter of 
8.6% is more than twice that of the GO case. However, 
the scatter in this model is mass-dependent; splitting the 
sample at 3 x 10^** M0, the low-mass end has ~ 13% 
scatter while it is only ~ 7% on the high-mass end. 

The opposing statistical shifts in /icm and effec- 
tively cancel when combined to form the thermal SZ 
signal. Not only does the mean Y — M relation in 
Figure [5] agree well between the two simulations, but 
both simulations have a similarly low scatter, ~ 12%. 
The small intrinsic sc atter is consistent w i th results from 
previ o us simulations (lEyrard et al. 119901: Ida S ilva et al] 



2000t iKav et al.l [2004 



Motletal 



20051 : iNaeai 2006 



Kravtsov et al T l2006l : lHallm~ et al As discussed 



below, the combination of low scatter in Y and the steep 
slope of the Y—M relation make the thermal SZ effect, or 
its X-ray equivalent, an excellent mass proxy for cluster 
surveys. Furthermore, as discussed above, the distribu- 
tion of deviations in ln{Y) at fixed mass is very close to 
a Gaussian distribution. 
In both models, the scatter in Tgi at fixed mass is higher 
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Fig. 10. — Graphical representation of the data comprising the covariance matrix for the two simulations at z = 0. Diagonal panels plot 
the distribution of deviations from mean mass scaling behavior, (s; — Si), where Si is the natural log of the i**^ signal. Shaded histograms 
show PH results and soUd lines show GO data. Each off-diagonal panel plots the normalized deviations, (s^ — Si)/cri, for an pair of 

properties. The lower triangle shows PH and the upper triangle shows GD behavior. 



than the scatter in at fixed mass. The difference is 
slight in PH, but over a factor of two in GO , a reflection of 
the larger amount of cool substructure in the latter treat- 
ment. The fact that Tgi is sensitive to the gas physics 
treatment means that opportunities exist for constrain- 
ing gas physics based on high-quality X-ray spectroscopy 
of large samples. Such a study could be provided by the 
proposed WFXT sateUite^. 

After NFW concentration, the X-ray luminosity has 
the highest level of scatter at fixed mass in both the 
PH {ainL = 0.19) and GO (aini = 0.28) simulations. 
These values are consistent with the core extr acted value, 
T = 0.27 ± 0.06, observed for REXCESS (|Pratt et all 
[2008h . 



4.3. Off-Diagonal Elements of the Correlation Matrix 
^ http://wfxt.pha.jhu.edu/ 



There is much information about the physical processes 
driving the ICM encoded in the correlation matrix of 
Table [B] For cosmological studies, it would be useful 
to identify pairs of cluster properties whose correlation 
coefficient is insensitive to gas physics modeling. High 
values of the correlation coefficient may point to pairs of 
halo properties that evolve on similar time scales during 
merger events. Below in §sec:imp, we show that high 
signal correlations can improve halo mass selection by 
joint signals considerably. 

Overall, the signal pair correlations tend to be higher 
in GO than in PH , but detailed differences warrant closer 
inspection. Considering the virial scaling as an anchor, 
we begin with consideration of correlations with ctdm- 
We then examine SZ and X-ray observables, and point 
out extreme values for both cases. 

The covariance between ctdm and Tm at fixed mass is 
an indicator of halo virialization. At redshift zero, this 
correlation is much higher, C = 0.56 in GO than in PH , 
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TABLE 5 

Scatter and Distribution of Deviates at Redshift Zero^ 



Q- 1 

bignai 


Scatter 


2.27570 


15.86570 


50% 


84.13570 


97.72570 


Gctussi&n 




-2.00 


-1.00 


0.00 


1.00 


2.00 


rn 




-1.66 


-0.90 


-0.11 


0.85 


2.43 


7ICM 




-2.14 


-0.93 


0.03 


0.93 


1.93 


PH Tm 


0.058 ± 0.002 


-1.90 


-0.84 


-0.04 


0.75 


2.11 


PH Tal 


0.069 ± 0.001 


-2.17 


-0.93 


0.04 


0.90 


1.88 


PH y 


0.125 ±0.002 


-2.04 


-0.97 


0.04 


0.90 


1.95 


PH Lbol 


0.193 ±0.002 


-1.98 


-0.98 


0.004 


1.00 


1.94 


PH Q 


0.116 ±0.001 


-1.89 


-0.99 


-0.04 


0.97 


2.19 


PH c 


0.300 ± 0.008 


-1.27 


0.43 


0.13 


0.89 


1.31 


GD (TDM 


0.042 ± 0.001 


-1.65 


-0.88 


-0.12 


0.77 


2.54 


SO /iCM 


0.036 ± 0.001 


-2.12 


-0.98 


0.03 


0.99 


1.90 


GO Tm 


0.102 ±0.001 


-2.23 


-1.00 


0.12 


0.93 


1.69 


GO T,i 


0.219 ±0.002 


-2.32 


-1.06 


0.19 


0.96 


1.49 


GO y 


0.123 ±0.001 


-2.17 


-1.00 


0.10 


0.95 


1.75 


GO Lbol 


0.282 ±0.003 


-2.19 


-1.03 


0.10 


0.97 


1.70 


GD Q 


0.109 ±0.001 


-2.59 


-0.91 


0.19 


0.85 


1.55 


GO c 


0.280 ± 0.008 


-2.42 


-0.97 


0.13 


0.86 


1.28 



" The distribution of deviates shows locations, in terms of normalized deviates, 
(5/(7, at which the ranked distributions reaches the listed percentiles. Values 
for a Gaussian distribution are listed in the first row. for each signal. Error 
estimates on the scatter come from bootstrap resampling. 



TABLE 6 

Correlation Coefficients at Redshift Zero ^ 



Signal 


0"DM 






/iCM 


y 


-^^bol 


Q 


c 


0"DM 




0.55 


0.81 


0.28 


0.54 


0.51 


0.17 


0.19 




0.35 




0.85 


0.48 


0.97 


0.67 


0.38 


0.49 


Ta 


0.86 


0.5 




0.42 


0.83 


0.67 


0.47 


0.64 


/iCM 


-0.10 


0.42 


0.37 




0.69 


0.60 


0.32 


0.37 


y 


0.079 


0.74 


0.62 


0.88 




0.73 


0.40 


0.51 


-^^bol 


0.26 


0.50 


0.73 


0.76 


0.78 




0.65 


0.70 


Q 


0.32 


0.029 


0.56 


0.15 


0.12 


0.59 




0.71 


c 


0.15 


0.053 


0.39 


0.29 


0.26 


0.51 


0.64 





The redshift zero correlation coefficients, with the results from 
the PH simulation in the lower triangl e a nd the results from the GO 
simulation in the upper, as in Figure [Tol Uncertainties from boot- 
strapping resampling are on the order of 0.01 and are not shown. 



where it is 0.35. The hales in the GO simulation are gov- 
erned by gravitational effects only, so and ctdm excur- 
sions track each other closely. In PH , however, the pre- 
heating raises the sound speed in the halos, making the 
thermalization due to mergers less pronounced, thereby 
diminishing (though not eliminating) the coupling of Tm 
and (TDM deviations. 

The ICM mass fraction behavior in the two simula- 
tions is quite different, as noted above, and this differ- 
ence is apparent in the covariance of ctdm and /icm- In 
the PH simulation, the two properties are anti-correlated, 
C ~ —0.10 at redshift zero. In fact, this is the only 
negatively correlated pair of signal deviations exhibited 
by the models. It is likely that this negative correla- 
tion is driven by the behavior of mergers. Halos in the 
early stages of a merger will have a higher ctdm at fixed 
mass than the mean relation. If, during these merger 
events, the collisionless dark matter accretes faster than 
the baryons (since the extended PH gas envelopes of ac- 
creting satellites will be more easily ablated during the 
merger encounter) , then the ICM mass fraction will be lo- 
cally depressed, driving an anti-correlation between ctdm 
and /icm- We note that these effects are absent in the 
GD case, since the correlation of that model is positive, 
C = 0.28. 



The concentration and ICM emission measure are 
measures that are sensitive to form ation history 
(|Wechsler et al.l 120021 : iBusha et al.ll2007f ) and substruc- 
ture driven by merging. The correlations between these 
parameters and (Tdm are surprisingly modest, ~ 0.2, per- 
haps indicating that formation history is a more impor- 
tant driver, compared to recent mergers, for these mea- 
sures. There is not a strong sensitivity to gas physics in 
the c— Ctdm correlation, but the PH Q — ctdm correlation of 
0.32 is significantly larger than the GO value. In fact, the 
only two signals with lower ctdm correlations in GD than 
in PH are Tsi and Q, both measures that are sensitive to 
small regions of cool, dense gas. 

Both Y and Lbol have a weaker correlation with ctdm 
in PH than in GO . As 1" cx /icMTm, we expect that we 
can describe the scatter in Y at fixed mass as 

ay — cr'^ + CTf + 2CTfcrT<Jf. (15) 

Examining Table EJ we see that, because the scatter in 
Y is largely dominated by T,„ in the GO model and /icm 
in the PH case, the value of the correlation coefficient, 
Ct/, is not an important factor. Our value for the scat- 
te r in Y is higher than that measured in the simulations 
of iKravtsov et all (f2006.) . as they see an anti-correlation 
in the X-ray inferred deviations of gas mass and temper- 
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ature at fixed mass. In the intrinsic measures of Table [6l 
we find positive correlation between /icm and Tm with 
values of 0.42 (PH ) and 0.48 (GO). 

It is worth noting the interesting case of C = 0.08 be- 
tween Y and ctdm in the PH simulation. Algebraically, 
this is unsurprising due to the positive correlation be- 
tween ctdm and Tm and the negative correlation with 
/icm- This lack of correlation suggests that SZ surveys 
should produce cluster samples that unbiased with re- 
spect to dynamical state. Similarly, the low correlation 
between Y and Q, C = 0.12, shows that SZ surveys 
should not be strongly biased by gas dumpiness. 

When metals are ignored, the bolometric luminosity 
scahng, L oc ficM^QT^^^, implies that the scatter in 
luminosity at fixed mass follows 



cr^ = 4ct^ +aQ 



4* 



fQ 



2* 



fT 



QT- 



(16) 



For the PH case, the scatter in ICM mass fraction is the 
primary contributor, responsible for 90% of the variation 
in Lboi- In contrast, /icm variations account for only 
7% of the scatter in the GO case, where the majority 
contributors are variations in Tsi and Q. Our r esults are 
consistent with the work of iBalogh et al.l (|2006f ). who use 
an analytic model to show that variation in halo structure 
cannot completely account for the observed variance in 
Lboi at fixed mass. 

In both simulations, the correlation between Tm and 
Ts\ at fixed mass is very high. Although a given halo 
will not have the same T,„ and Tgi, the high correlation 
coefficients indicate that the two temperature measures 
similarly trace the thermal state of a halo. The correla- 
tion is higher in the PH simulation than in the GO , and 
cool cores in the latter model also drive higher correla- 
tions between Tsi and Q than between Tm and Q. The 
correlations between /icm and temperature measures are 
nearly constant between the two simulations, with values 
in the range 0.4 — 5 for both Tm and Tgi- 

The highest measures of correlation are between Tm 
and Y in the GO case (0.97), and between /icm and 

Y in the PH run (0.88). The SZ-X-ray correlation of 

Y and £boi is also large, ~ 0.7, in both treatments. 
This robust behavior is promising for cr oss-calibrations 
of future, combined X -ray and SZ survevs lVounger et al.l 
(j200l);[Cunh3 dOOl). 

Both simulations also have a significantly non-zero cor- 
relation coefficient between Lboi and Tgi, 0.73 in PH and 
0.67 in GO . As shown bv Nord et all (2008). the X - T 
scaling relation expected from X-ray flux-limited surveys 
is sensitive to the value of Ctl, and studies of current 
and future samples will be able to place limits on this cor- 
relation, us i ng tec hniques similar to those employed by 
iRozo et all ()2009D to constrain the correlation of mass 
and X-ray luminosity at fixed optical richness for the 
SDSS maxbcg sample. 

4.4. Redshift Evolution of the Signal Covariance 

We plot the time evolution of the signal scatter in Fig- 
ure lll[ with error bars come from bootstrap resampling. 
Going back to z = 2, few halo properties show any evo- 
lution with redshift. In the PH case, the scatter in baryon 
fraction slightly increases near z — 2, causing an increase 
in the scatter of Y and L. On the other hand, the scatter 
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Fig. 11. — Evolution of the scatter in eight bulk cluster proper- 
ties. 



0.5 




0.5 - 



Fig. 12. — Four pairs of signals which show little to no evolution 
in the correlation coefHcient. The black, sohd line denotes the PH 
simulation, and the red, open points the GO simulation. 

in emission measure decreases at higher redshift. Note 
that, at redshift two, there are only 62 halos in the PH 
simulation above our mass cut of 5x 10^^h~^MQ. At red- 
shift zero, there is a higher scatter at the low-mass end 
of our mass range, as seen in Figure [2l Since the halos at 
redshift z = 2 are only slightly over the mass cut, their 
higher scatter may reflect this mass dependence rather 
than pure redshift evolution. 

Moving on to the off-diagonal elements, most pairs of 
signals show little evolution in the correlation coefficient 
with redshift. Figure [12] shows four typical pairs with 
little evolution. Even as the physical density of halos 
changes with redshift, we see that the interplay of Y and 
L or Tsi and /icm does not change. We see particularly 
little evolution in pairs of signals in the PH simulation. 

We do see evolution in a few signal pairs, notably be- 
tween (TDM and other signals in the GO simulation. Sev- 
eral pairs of signals - such as ctdm — Tsi and /icm — Y- 
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Fig. 13. — Four pairs of signals which show some degree of evolu- 
tion in the correlation coefficient, particularly in the GO simulation. 
The blaclc, solid line denotes the PH simulation, and the red, open 
points the GO simulation. 

evolve in the GO simulation, but not in the PH simulation. 
Furthermore, for ctdm — ^bot, we see evolution in both 
simulations, but in opposite directions. The evolution of 
these pairs is shown in Figure [T51 

The correlation between ctdm and the SZ and X-ray 
signals increases towards high redshift in the GO simula- 
tion. We speculate that this is due to the increase in Q 
with time at fixed mass in the GO simulation, as shown 
in Figure [S] The gas in the GO simulation develops more 
substructure with time, and this tends to decouple the 
SZ and X-ray measurements from the dark matter ve- 
locity dispersion, (Tdm- Although Q also increases with 
time in the PH simulation, its normalization at all red- 
shifts is much lower than in the GO simulation. Hence, 
the substructure of the gas does not contribute much to 
the Tsi or L measures in the PH simulation, and there is 
little evolution of their correlation with (tdm- 

4.5. Implications for Multi-Signal Mass Selection 

At fixed signal, the scatter in halo mass is cr^i = Ci/ai, 
where i labels the particular signal and at is the slope of 
that signal-mass relation. From the analysis above, we 
compute the mass scatter for a subset of z = signals 
and present the data in the first two columns of Table 
[7l We see that Y provides the best mass selection un- 
der both physical treatments, with scatter ~ 7%. For 
the PH case, the X-ray luminosity is quite good, with a 
10% scatter, while /icm is the worst selector, with the 
weak /icm — M slope producing a 28% scatter in mass. 
We will see below that this large scatter can actually be 
used to improve selection when paired with more precise, 
correlated signals. 

When selecting halos using multiple signals, the mass 
variance is = (a^^'~^a)~^, where a is a vector of 
slopes with elements a^. We consider the two-signal case, 
and let r = C12 be the (traditional) correlation coefficient 
between signals 1 and 2. 

Given a pair of signal measurements, the mass selection 
is log-normal with variance 

= (1 - r'){a-J + a-J 2ra;la-J)-\ (17) 
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Fig. 14. — For a range of correlation coefficients, we plot the 
dispersion in mass achieved by combining a pair of signals (1 and 
2) against the ratio of their dispersions. Different lines vary the 
correlation coefficient r between the pair of signals at fixed mass. 
Given a measurement of the first signal, a second measurement 
always improves the mass scatter, except in the degenerate case, 
""mz/i^'/jI = 1 S'lid r = \. While negatively correlated signals yield 
more improved mass selection than positively correlated pairs, the 
improvement becomes independent of sign when the second signal 
is much noisier than the first, (T^2/o"(il 3> 1. 



The improvement in mass selection to be gained by a 
pair of signals, relative to a single measurement, is dis- 
played in Figure [m Here, we plot S/cr^i as a function 
of the scatter ratio (0-^2 /c^i) for several values of the 
correlation coefficient. If the second signal is a better 
mass proxy, a^a <C ct^i, then it dominates the selec- 
tion. Combining signals with comparable mass selection, 
<^iii ^ cT/jIi can result in anything from no improvement 
in the degenerate case (r ^ 1) to the \f2 improvement 
for uncorrelated signals (r = 0) to dramatic improvement 
in the anti-correlated case (r —* —1). Surprisingly, when 
an intrinsically "noisy" mass proxy is added, a^a S> cr^i, 
one still achieves significant improvement in mass selec- 
tion as long as the signal pair correlation is large in an 
absolute sense. 

We evaluate equation p7|) for a set of observable signal 
pairs and present the resultant mass scatter values in the 
right-hand columns of Table [7l In the PH simulation, 
the correlation coefficient between L and y is r = 0.78, 
and the scatter in mass for the pair of signals is 6.9%, a 
mild improvement over the scatter using only Lboi (11%) 
or only Y (7.5%). Given the high covariance between 
Lboi and Y ^ future multiwaveband surveys that join SZ 
and X-ray detections should have a strong mass selection 
properties. 

Finally, we note that the best mass selection comes in 
the PH case from combining Y with /icm- The combina- 
tion of a strong correlation, r — 0.88, and, especially, the 
relatively large degree of scatter in mass at fixed /icm 
(0.28) produce a scatter in mass of only 4.1% when mea- 
surements of these two signals are joined. We caution 
that this analysis does not take measurement errors into 
account, and it is clear that high data quality will be 
important to realize this level of mass selection. 
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TABLE 7 

Mass Scatter at Redshift Zero ^ 



Signal 


0"DM 




/iCM 


Y 




o"DM 


0.12\0.12 


0.12 


0.12 


0.075 


0.12 


Tel 


0.10 


012\0.38 


0.35 


0.050 


0.26 


/iCM 


0.11 


0.12 


0.28\0.12 


0.054 


0.21 


Y 


0.062 


0.069 


0.041 


0.069\0.075 


0.066 




0.09 


0.10 


0.09 


0.069 


0.10\0.26 



^ The diagonal elements give the mass scatter for individual signals for 
the PH\G0 cases. Off-diagonal elements give the halo mass scatter for 
signal pairs, with the P H c ase in the lower triangle and the GO case in 
the upper, as in Figure [Tol . 



5. CONCLUSION 

We analyze scaling relations for multiple properties of 
massive halos taken from a pair of gas dynamic simula- 
tions with different gas physics treatments. Our sam- 
ples contain tens of thousands of halos with masses 
M200 > 5 X 10" h-^ Mq at redshifts z < 2. The phys- 
ical treatments of gravity only (GO) or preheating (PH) 
are both highly idealized, but we show that the latter re- 
produces the scaling relation behavior of core-extracted 
X-ray measures of local clusters. 

The dark matter velocity dispersion scales with mass 
and redshift according to self-similar expectations, indi- 
cating that the virial theorem is respected regardless of 
gas treatment. However, the gas behavior in both treat- 
ments differs from self-similarity. The deviations in the 
GO case tend to be small and are related to the mass- 
limited sample definition. 

The entropy injection of the PH model drives more sub- 
stantial deviations from self-similar scaling. At z = 0, 
the ICM mass fraction varies with mass in a manner 
roughly consistent with observed measurements of local 
clusters. We find that /icm requires a quadratic fit in 
log-mass, and mild curvature in the logarithmic scalings 
of Y and Lboi as a function of mass is also evident. While 
the ICM in PH halos is lower in mass compared to the GO 
case, it is also slightly hotter. The effects on /icm and 
Tm nearly cancel when combined to form the thermal SZ 
signal, leading to nearly identical Y-AI scaling relations 
above 3 x 10^'* for both PH and GO cases at low redshift. 

The ICM mass fraction at fixed mass declines weakly 
with redshift, by 10% in 5 x lO" /i'^ Mq halos at z = 1 
and with larger reductions at lower masses. While Chan- 
dra observations of optically-selected clusters in the RCS 
survey show evide nce of reduced gas fractions at z — 1 
(jHicks et al.ll2006l ) , further work is needed to address the 
quantitive level of agreement between the observations 
and PH model expectations. The PH baryon fraction evo- 
lution drives departures from self-similar predictions in 
the Y-M and L^oi-M relations; slopes tend to steepen 
slightly and the normalization at 10"'^^ Mq is lower 
than self-similar expectations at high redshift. 

We present the first systematic investigation of prop- 
erty covariance in the computational samples of massive 
halos. The data generally support a multivariate log- 
normal form for the joint distribution of signals at fixed 
mass and redshift. All measures depart somewhat from 
an exact gaussian form, but the deviations in gas mea- 
sures are smaller for the PH model due to the suppression 
of substructure caused by the preheating. Most signal 
pairs exhibit positive correlations, with the lone excep- 



tion of —0.1 between (Tdm and /icm in the PH case. The 
thermal SZ signal displays a robust 13% scatter that is 
strongly correlated with variations in both ICM gas mass 
and temperature, with /icm dominating in the PH case 
and Tm being more important in the GO treatment. 

Combining multiwavelength observations offers an op- 
portunity to improve selection of clusters by their intrin- 
sic mass. We derive the mass variance of signal pairs and 
show that combining strongly correlated signals always 
improves mass selection, even when one of the signals by 
itself is a comparatively poor mass proxy. The combina- 
tion of thermal SZ and ICM mass fraction in the PH case 
selects halo mass with just 4% intrinsic rms scatter. 

Identifying the root causes behind the terms in 
the covariance matrix is a consider able task that 

we le a ve for future wo rk. Mergers (jRoettiger et al.l 

11997"; 'Ricker & Sarazin' '2001|) and assembly bias 
(Boylan-Kolchin et al. (2009) and references therein) will 
surely play important roles for many cluster signals. 
Studies of merger history behavior will shed light on sur- 
vey selection properties, particularly potential biases re- 
lated to the dynamical state of a halo. 

The simple treatment of baryon physics in our simu- 
lations limits our investigation to the hot, thermal ICM 
of clusters. More extensive physical treatments that in- 
corporate galaxy and supermassive black hole formation 
and other physics such as MHD and non-thermal plasmas 
will ultimately extend the set of observable halo signals 
into the optical/NIR and radio wavebands. 
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